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I. INTRODUCTION 

In conventional quantum mechanics, one demands that 
the Hamiltonian H generating the time evolution has a 
real spectrum and that the corresponding time evolution 
operator U is unitary. These conditions are fulfilled if the 
Hamiltonian is Hermitian, i.e. H = H\ which is usually 
considered an axiom of quantum mechanics. However, 
the condition of Hermiticity can be weakened. In the 
class of so-called PT-symmetric Hamiltonians, one can 
ensure real eigenvalues and a unitary time evolution even 
for explicitly non-Hcrmitian Hamiltonians [2, 3]. A thor- 
ough review of the foundations of PT-symmetric quan- 
tum mechanics can be found in [1]. 

In the following we will investigate entanglement phe- 
nomena in bipartite systems within this framework. An 
entangled state is a quantum state where two or more 
degrees of freedom are intertwined, so that they are not 
independent anymore. In this context a couple of histori- 
cal discussions took place and gave deep insights into the 
nature of quantum mechanics, like the Einstein-Podolsky- 
Rosen paradox [8]. These states have a wide range of 
applications, for example in quantum information theory 
and quantum computing. The question of entanglement 
generation and entanglement efficiencies for conventional 
quantum mechanics was addressed earlier in [7]. 

Using PT-symmetric quantum mechanics to describe 
entanglement phenomena was also done in [13]. Espe- 
cially for the case of bipartite systems, relations for the 
degree of entanglement of given states and the entangle- 
ment capability of certain systems could be found. Al- 
though we confirmed many of these results, we have some 
different and new findings. For a particular initial state, 
as well as for general states, we give relations between 
the efficiency of the system to generate entangled states 



and the parameters of the Hamiltonian describing the 
dynamics of the system. 

We will first introduce a measure for entanglement and 
a certain class of PT-symmetric Hamiltonians. We then 
quantify the degree of entanglement of an arbitrary and 
of a generalized Einstein-Podolsky-Rosen state. Subse- 
quently, we are dealing with the question of entangle- 
ment generation of a particular PT-symmetric state and 
generalize for arbitrary states. The question we try to 
answer is how to characterize the rate of entanglement 
generation for this class of systems. 



II. BIPARTITE SYSTEMS 
A. A measure for entanglement 

Let ~Hi=i,...,N denote a set of Hilbert spaces. We call a 
state of a composite system H = rl\®- ■ -^Hn entangled, 
if there is no decomposition of the form |W) = |xi)<8>- • -<8> 
\Xn) with suitable \x%) £ In the following we restrict 
ourselves to bipartite systems, i.e. N = 2. 

A measure of entanglement are the entropies 

E(V) = - trj ( Pl log 2 px) = - tr 2 (p 2 log 2 p 2 ) , (1) 

where p\ — tr 2 p and p 2 — tri p are the reduced den- 
sity matrices, p = \^f) (^f\ is the density matrix itself and 
trj denotes the partial trace over the i th subsystem [12]. 
Here E($) g [0, 1] and E($) = if and only if |W) € U 
is not entangled. 



B. General entanglement content 



Consider the Hamiltonian 



H 
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with s 2 ^ r 2 sin 2 0. Observe that H is in general not 
represented by a Hermitian matrix, W 7^ H. But it 
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obeys PT-symmetry with 



V = 



1 

1 



(3) 



and T complex conjugation [4] (see also erratum [5]), 
i.e.[H,VT]=Q. 

Define sint^ = r/s ■ sin0 with ip S [—tt/2, tt/2]. The 
simultaneous eigenstates of H and VT are given by 



-) 



1 



yJ2 cos ip 



y/2 COS<^ 



e -iW2 



(4) 



with eigenvalues E± — r cos ± s cos y>. These eigen- 
states are orthonormal with respect to the positive CVT- 
inner product (^f\^)cvT = Wcpt ' |$)> where (^\cvt = 
(CPT|^)) T , see [4]. For this Hamiltonian, the C operator 
has the form 



C = 



i tan if sec p 
seep — i tan (/? 



(5) 



with [iJ,C] = 0. 

In the following we will consider the composite system 
% = Hi x H 2 °f Hilbert spaces Hi = i,2 with dynamics 
governed by the Hamiltonian H = Hi £g> H 2 , where Hi 
is of the form of (2) with r i ,s l ,Q i S R, sf ^ r 2 sin 2 6i, 
sin i/3j = ri/sj-sin 6; and \4>±)i G Hi being the eigenstates 
of i/i for i = 1,2. Our respective C operator reads Ci (g) 
C2, where is of the form of (5) with corresponding 
parameters. As 



[Hi ®H 2 ,VT® VT] = [Hi ®H 2 ,Ci® Ca] 
= [Ci<g>e 2 ,PT<g>P7] 
= 



(G) 



holds, H obeys VT symmetry and we can define a CVT- 
inner product as above. The decomposition of the op- 
erators means the same direct-product form of the met- 
ric operator 77 in [14]. According to proposition 2 and 
proposition 3 in the same paper, this permits a proper 
quantum mechanical description of the bipartite system 
with a unitary time evolution. 

Now consider the density matrix p = \^f(i))(^f(t)\cvT 
of a system in the state 

|*(t)) = a{t)\<t>+)i ® \4>+) 2 + P{t)\<t>+)i <8> |0-)a 

+ 7 (t)k-h <8> l^+>2 + ® |0-) 2 (7) 

with Y$>(t)) € The eigenvalues of pi(t) = tr 2 p(t) are 
given by A±(t) = ^iV 1 ~ S W withS(t) = 4|a(*)<f(*)- 
/3(t)7(t)| 2 , which is a simplification of the result in [13]. 
Hence the entanglement content is 

E(t) = E(*(t)) 

= -X+(t) log 2 A+ (*) - A_ (t) log 2 A_ (t) . (8) 

Note that \^f(t)) only separates if a(to)6(to) = f3(t )j(to) 
for t = to. 




Figure 1: Entanglement content E^ ) as a function of <pi 
and Lp2- 



C. The Einstein- Podolsky- Rosen state 

Consider the Einstein-Podolsky-Rosen state from con- 
ventional quantum mechanics and normalize them with 
respect to the CPT-inner product: 



= Vcos Tp~i f J ' ^ = V cos <Pi ( 1 J 



Overall normalization yields 



|*-> = «(|t>l® lt>2) 



(9) 



(10) 



with k = [2 (1 — sin ipi sin <^> 2 )] 



-1/2 



The eigenvalues of 



the reduced density matrix of the first subsystem are 
given by 



1 sin (pi — sin ip 2 

2 2 (1 — sin ipi sin <y5 2 ) 



(11) 



Our results are in disagreement with [13] mainly due to 
the author's use of non-CPT-normalized states. 
We can consider the entanglement content 



£(* ) = -A+ log 2 A+ - A_ log 2 A_ 



(12) 



as a function of ipi and tp 2 , i.e. of the Hamiltonians Hi 
and H 2 . The result can be seen in figure 1. Note the 
cases <pi — (p 2 or ipi + ip 2 = ±tt, where (10) has the form 
of a 'PT-symmetric Bell state. 
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Figure 2: Entanglement content E(t) (solid) and entangle- 
ment rate r(t)/r~ 1 (dashed) of |\£(t)) with r = 2ft/ (uiiuz). 



III. ENTANGLEMENT GENERATION 

A. Entanglement capability 

We investigate the question how to increase the capa- 
bility of a system to entangle states. More precisely we 
want to understand the dependencies of the entanglement 
capability T(t) = dE(t)/dt of the parameters of the sys- 
tem. In order to determine the time evolution operator 
of the Hamiltonian Hi (g> H2 (see section II B) define 



rii = — (sj, 0, i^sinOi) 1 



2si cos tpi 



(13) 



Rewrite H as 



H = (r t cos Oil + y nt ■ 



r 2 cos6 2 l+y^-^) , (14) 

with a = {0 1 , (J2, (^3)^ 3 <Tj denoting the Pauli matrices 
and 1 the 2x2 identity matrix. Expanding the expression 
for H yields four terms of which only one can generate 
entanglement. We restrict ourselves to this term and 
define H = W1CJ2/4 • (rii ■ a ® ■ a) resulting in a time 
evolution operator U(t) — exp(—\Ht/h) given by 

f/(t)=cos(^) 

-isin(^)(^.^®^.^). (15) 

Consider now a simple initial state, namely \*5f(t — 
0)> = If)® It), and apply U(t) to find 

|¥(t))=a(t) |t)®|t} + 0(t) It)® It) 

+ 7(*) It)® |t> + *(«) It)® It) (16) 

with 

/ uiUJ2t\ 4irir 2 sin 0i sin 02 /wit^iA 

a{t) = cos \w ) + — ^ sin ) ' 



4s 2 rism6i . ( wiu$t 

p{t) = sin — 

W1W2 V 4ft 



7(*) 



4siT2sin02 ( UiWit 
= sm ~n~ 

LO 'iU>2 V 4fi 



4isis 2 . / 

c5(t) = sm — — 

a;ia;2 \ 4ft / 



(17) 



Up to here our results are in agreement with [13]. We 
now find the entanglement content to be given by 

E(t) = -\ + {t) log'2 A+(t) - \-{t) log 2 A_(t) (18) 



with 

A ± (t)= 2±2 cos(^ — j 

We find the entanglement rate to be 
dE(t) 



(19) 



T(t) = 



dt 



= — ^ — jlog 16 cot (20) 

The maximal entanglement capability is r max = 
max t r(£) = 0.4781 Wia;2 /ft- Therefore by changing u>i 
due to an adjustment of the parameters of the Hamil- 
tonians Hi, we can control r max . The typical time de- 
pendency of the entanglement content and entanglement 
rate can be seen in figure 2. 



B. Efficiency of general systems 

If we have a given Hamiltonian one may ask how to 
maximize the entanglement rate of the system. We want 
to generalize some results of conventional quantum me- 
chanics addressed in [7] to the PT-symmetric case. 

Consider an arbitrary PT-symmetric two qubit sys- 
tem. Using the (Sc/imidf-dccomposition theorem we 
rewrite an arbitrary state \^(t)) <G Hi x H2 as 

!*(*)> = Vp$)\vt) ® \xt) 

+ e' a ^p{f)U)^\xt) (21) 

with p£ [0, 1] as ^c/imidi-coefficient and (tpt\<Pt )cpt — 
(Xt\Xt)cvT — as Schmidt-vectors. The entanglement 
content is given by 

E(t) = -p(t) log 2 p(t) - (1 - p(t)) log 2 (1 - p(t)) (22) 

and the entanglement rate factorizes in two terms, 
i.e. T(t) = dE(t)/dp(t) x dp(i)/dt, where 



dE(t) 2 
dp(t) ~ log2 



arctanh (1 — 2p(t)) . 



(23) 



After choosing the phase a appropriately the evolution of 
the Schmidt-coefUcierit is determined by the differential 
equation 
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dp(t) 



VK*)(i-f(*)) 



dt ft 

x | {(<Pt\cvT ® (x*|ct>t) (bt X ) ® |Xt")) 



(24) 



This relation is known from conventional quantum me- 
chanic, but also holds for "PT-symmetric systems. Define 

hi! = max | ((<p\cvT ® (Xlcvr) H (|^> ® | X X )) I- (25) 
IMI =1 

llxll=i 

Then the time evolution of p(f ) for an optimally prepared 
setup, i.e. a setup forcing the qubit states to be optimal 
at every instant of time 



I {(<Pt\cPT ® (Xtlcvr) h (lv»*-> ® IXt")) I - 

for all t, follows from the differential equation 
dpopt(t) 



(26) 



dt 



2fty^ op t(t)(l-popt(t)). (27) 



We find Popt(t) = sin 2 (fit + <5o) with an integration con- 
stant Sq , which is in agreement with results from conven- 
tional quantum mechanics in [7] (see also [13]). Hence we 
can characterize the entanglement rate of an optimally 
prepared setup completely in terms of f2 via 



r op t(t) 



d P (t) 



dpopt(t) dE(t) 
dt 

2Q 
log2 



(28) 



Popt(t) 

arctanh (cos (2fit + 2S )) sin (2fit + 2S ) 



and find for the maximal rate r max = max t r op t(i) = 
1.9123 f2. The parameter is the critical value one needs 
to maximize to efficiently entangle states. We remark 
that this result also holds in conventional quantum me- 
chanics, where is defined with conventional conjugates 
instead of CPT-conjugates. 



IV. CONCLUSIONS 

For the case of a state from a bipartite system we de- 
termined the degree of entanglement and saw the emer- 
gence of symmetrical patterns (see figure 1) in the case 
of the Einstein-Podolsky-Rosen state. We quantified the 
capability of a given ■pT-symmetric system to generate 
entangled states in terms of the parameters of the Hamil- 
tonian. Their ability to entangle states can be described 
by the parameters CJj=i 2 in (13) or in general by f2 in 
(25). Many relations are similar to the corresponding 
ones in conventional quantum mechanics after replac- 
ing the usual inner product with the CP7~-inner prod- 
uct. However, these results are not obvious and need 
to be checked. For example, the recent discussion of 
the quantum brachistochrone problem showed that VT- 
symmetric quantum mechanics can give some surprising 
results [6, 9-11]. 
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